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Abstract 

In this paper we investigate the boundary value problem 

div(7Vii) = in fi 

u — f on dSl, 

where 7 is a complex valued L°° coefficient, satisfying a strong ellipticity 
condition. In Electrical Impedance Tomography, 7 represents the admit- 
tance of a conducting body. An interesting issue is the one of determining 
7 uniquely and in a stable way from the knowledge of the Dirichlet-to- 
Neumann map A T . Under the above general assumptions this problem is 
an open issue. 

In this paper we prove that, if we assume a priori that 7 is piecewise 
constant with a bounded known number of unknown values, then Lipschitz 
continuity of 7 from A 7 holds. 

1 Introduction 

In this paper we investigate stability for the inverse problem of electrical impe- 
dance tomography. More precisely we consider the following problem: let u G 
H x (£l) be the solution to 

div (7VU) = in ft . , 

u = f on <9f2, ' 

where fl C M", n > 2 is a bounded connected domain, 7 is a complex valued 
function representing the admittivity coefficient, it is bounded and satisfies the 
ellipticity condition 5R7 > A -1 > a.e. in Vl and / G H l l 2 (dVl). 
The Dirichlet to Neumann map A 7 is the operator A 7 : H l / 2 {d£l) ->• H^^idfl) 
given by 

ov\ aa 
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where v is the exterior unit normal vector to dQ. 

The mathematical formulation of the inverse problem of impedance tomog- 
raphy is to determine the admittivity 7 from the knowledge of the Dirichlet to 
Neumann map A 7 . 

This problem has several important applications in fields like medical imag- 
ing and nondestructive testing of materials. We refer to the review papers by 
Borcea ( |Bo] ) and to f [CIN] ) for a wide bibliography on relevant examples of ap- 
plications. We want to point out that equation JT]) also appears in the study of a 
model for electrical conduction in biological tissues as the asymptotic limit of an 
elliptic equation with memory when subjected to periodic Dirichlet boundary 
conditions (see [AABGj and [AABG2] ). 

For n > 3 the uniqueness result by Sylvester and Uhlmann ( |SUj ) obtained 
for real conductivities applies also to the complex case (cfr. [Bo]). For n = 2 
the first contribution on the unique determination of 7 by A 7 was given in [F] 
where the author proved uniqueness provided the imaginary part of 7 is suffi- 
ciently small. In 2008 Bukhgeim in [Buj generalized this result to an arbitrary 
sufficiently smooth admittance. 

The problem of determining uniquely an arbitrary L°° admittivity from the 
Dirichlet to Neumann map is completely open even in the real case when n > 
3. The only result known is the one of Astala and Paivarinta who proved 
uniqueness of real L°° conductivities from the Dirichlet to Neumann map in the 
two dimensional case, (cf. [AP]). 

The main topic of our paper is to investigate continuous dependence of 7 on 
A 7 when the admittivity is an L°° function of a particular form. 

In general for arbitrary conductivities it is well known that this problem 
is severely ill-posed. If 7 is a real valued coefficient satisfying suitable a-priori 
smoothness assumptions, Alessandrini proved in [X] a log-type stability estimate 
for n > 3; the same type of stability was proved in |BFRj for n = 2 for Holder 
continuous conductivities. Such estimates are optimal (see [M])- Clearly one 
expects the same kind of ill-posedness also in the complex case. 

On the other hand, in many applications one has to disposal additional a 
priori information on the unknown function the might lead to better stability 
bounds. In |AV] Alessandrini and Vessella assume that the real conductivity 7 
is of the form 



where Dj are known disjoint Lipschitz domain and 7^ are unknown real numbers. 

Assuming ellipticity and C 1,Q regularity at the interfaces joining contiguous 
domains Dj and at <90 they prove Lipschitz continuous dependence of 7 on A 7 . 
The key ingredients in their proof are, on one hand the use of the Green's func- 
tion and its asymptotic behaviour near the regular interfaces, on the other hand 
the use of global C a regularity estimates of solutions and local C 1,a regularity 
estimates in a neighborhood of the smooth interfaces. 

In this paper we generalize the result in [AVj to the complex equation (fTj). 
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More precisely we show that if and are of the form 

AT 

7 (fc) (*)=E^ fc V(z), k = l,2 

with > A" 1 > for k = 1,2 and assuming that the interfaces joining 

contiguous domains contain a flat portion then 

Il7 (1) -7 (2) IU~(f2) < C||A 7 (d - A 7 (2)|| £ j- ffl/2( g 0))ff _ 1/2(an) -j, 

where C depends on fi, A and N and diverge to +00 exponentially as N — > +00. 

Our approach follows the one of Alessandrini and Vessella of constructing 
singular solutions and of studying their asymptotic behaviour when the singu- 
larity approaches the discontinuity interface. 

Observe that if 5R7 > A -1 > the complex equation ([I} is equivalent to a two 
by two strongly elliptic system with L°° coefficients. One relevant difference with 
the conductivity case treated in [AVj is that in the case of real L°° conductivities 
existence of the Green's function in f2 is guaranteed by the results contained 
in [LSW] while for equation (JXJ) , L°° admittivities and n > 3 the existence of 
the Green's function in the whole domain fi is not known due to the lack of 
a maximum principle and of De Giorgi-Nash type regularity estimates for this 
type of equations. 

We are able to bypass these difficulties observing that in order to derive our 
result it is enough to construct and to study the behaviour of singular solutions 
in a Lipschitz subset K (defined in Section 3) of a slightly enlarged domain Hq 
containing the smooth portion of the interfaces and determining its asymptotic 
behaviour near the interfaces. On the other hand in the domain K, using the 
estimates for elliptic systems obtained by Li and Nirenberg in |LN| . we have 
that solutions to equation ([1]) enjoy Lipschitz estimates and are C°° in each 
"strip" of K up to the flat interface. We present our analysis in the case n > 3 
although our result can be extended easily to the case n — 2. In fact the two 
dimensional case is in some sense easier to treat since in this case Dong and 
Kim in |DK] have proved existence, uniqueness and pointwise estimates of the 
Green's function in O. 

We want to point out that our result holds also if K contains less regular C 1,a 
interfaces (see Remark 5.1) but we think that the treatment of this case would 
only require tedious and long technicalities and calculations. On the other hand 
the flatness assumption of a portion of the interface is not too restrictive since it 
includes for example a partition of f2 with polyhedral domains Dj which appear 
in any numerical scheme used for the effective reconstruction of the admittivity. 
Besides, this assumption allows to derive Holder quantitative estimates of unique 
continuation of solutions to equation ([!} and consequently a better dependence 
of the constant C on N . 

The plan of the paper is the following: In Section 2 we introduce notation and 
the main assumptions and we state our main result (Theorem 2.1). In Section 3 
we collect all the results needed in order to prove Theorem 2.1. In Proposition 
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3.1 and Proposition 3.2 we state some known results concerning respectively the 
regularity of solutions of equation (J]) and the existence of the Green's function 
in the case of continuous admittivities. In the key Proposition 3.3 we prove 
the existence of singular solutions in K and we investigate their asymptotic 
behaviour near the flat discontinuity interface. In Theorem 3.4 we show that 
the solutions of ([1]), which due to the particular structure of 7 are piecewise 
analytic in il, can be extended analytically through the flat interfaces. This 
property allows us in Proposition 3.5 to derive optimal quantitative estimates of 
unique continuation for solutions of equation (fT|). In Section 4 we give the proof 
of Theorem 2.2. In Section 5 we give some final remarks about generalizations 
of our result and, finally, the Appendix contains the statement of Caccioppoli 
inequality (Proposition 6.1), the proof of Theorem 3.4, and the generalization 
of Alessandrini's identity to the complex case. 

2 Main result 

2.1 Notation and main assumptions 

For every x € R" let us set x = (x',x n ) where x' £ R 71-1 for n > 3. With 
Br{x) and B' R {x') we will denote respectively the open ball in E™ centered 
at x of radius R and the ball in M n_1 centered at x' of radius R\ Br(0) and 
B' R (0) will be denoted by Br and B' R . We will also use the following notations 
R\ = {{x',x n ) E R n : x n > 0}, Wl = {(x',x n ) £ I" : i„ < 0}, B+ = BrPiR^ 
and B R = B R n Rl . 

We will denote by D@ the derivative corresponding to a multiindex /3 = 
(/3i, ...,/?„) and by D^, the partial derivative corresponding to the multiindex 
j3' = (/?!,... , /3„_i, 0), while we will write for h € {1, ... , n}, for the partial 

derivative with respect to Xh and the partial derivative in direction v. 

Let f2 be a bounded domain in R™. We shall say that dfl is Lipschitz con- 
tinuous with constants tq, L > if VP € dfl there exists a rigid transformation 
of coordinates such that P = and 

ft n B ro = {(x', x n ) G B ro : x n > <j){x')}, 

where <p is a Lipschitz continuous function on B ro with 4>(0) — and 

||^Hc°.i(s; ) ^ Lr o- 

Our main assumptions are: 
(HI) £1 C is a bounded domain such that 

|n| < Ar%, 

and 

d£l is of Lipschitz class with constants rg, L. 
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(H2) The complex conductivity 7 satisfies 

^7>A~\ M<A in Q (3) 
for some A > 1, and is of the form 

N 

where 7,- are for j — 1, . . . , N unknown complex numbers and Dj are known 
open sets in R™ which satisfy the following conditions 

(H3) Dj, j = 1,...,N are connected and pairwise nonoverlapping such that 
Li^ =1 Dj = f), dDj, j = 1, . . . , N are of Lipschitz class with constants ro, L. 

We also assume that there exists one region, say Di such that dD\ n dVt con- 
tains an open flat portion Si. For every j £ {2, . . . , N} there exist ji, . . . , j^i G 
{1,...,N} such that 

D jl =D 1 , Dj M = Dj 



and, for every k = 1, . . . , M 



dD^ndDj, 



contains a flat portion S^ such that 

S fe c fl, Vfc = 2,...,M. 

Furthermore, there exists Pk G S^ and a rigid transformation of coordinates 
such that Pk = and 



Sfc n B ro/3 = {ie B ro/3 
Dj h nB ro/3 = {xeB ro/3 
D o k -i n B ro/3 = {xeB ro/3 



.r 



n = 0}, 

x n > 0}, 
x n < 0}. 



For simplicity we will call Dj 1 , . . . , Dj M a chain of domains connecting D\ 
to Dj. 

In the following we will introduce a number of constants that we will always 
denote by C . The values of this constants might differ from one line to the 
other. We will write explicitly which a priori parameters each constant depends 
on. 

Consider the problem 

div (7VU) = in 57 , . 

u = f on <9f2, 

where fl and 7 satisfy assumptions (H1)-(H3) and / G _ff 1 / 2 (<9i7). Observe 
that, by assumption ([3]), applying Lax-Milgram Theorem, there exists a unique 
solution u e H 1 ^) of problem (gj). 
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For / e H 1 ^ 2 (dfl) define the Dirichlet-to- Neumann map A 7 : H 1 / 2 ^) 
H-^ 2 (dn) given by 

du 

dv I an ' 



A 7 (/)=7- 



Let us notice that A 7 can be idcntihcd by the sesquilinear form on 7J 1 / 2 (<9i7) x 
H l ' 2 {dVt) defined by 

<A 7 /,V>= / 7Vm-VT7, V/,V G H^ 2 {dil) 
Jn 

where u is solution to problem Q and u is any function such that 

v \ a n = i>- 

Theorem 2.1 Let i7 satisfy assumption (HI). Let k = 1,2 6e fu>o com- 
ptec piecewise constant functions of the form 

N 

7 W (z) = E7?V(*)> 
i=i 

where 7^") satisfy for k — 1,2 assumption (H2) and -C^, j = 1,...,AT satisfy 
assumption (H3). 

T/iera t/iere exists a positive constant C — C[tq, L, A, n, N, A) smc/i i/iai 

|| 7 W - 7 ( 2 )|| iM(n) < C||Ai - A 2 || jC(jr i/» ( en),ff-i/»(an)). (5) 
where Aj = A„(&) /or fc = 1,2. 

3 Preliminary results 

Let us now state some results that will be useful in the proof of Theorem 12.11 

The first ingredient is a regularity estimate for solutions to the admittance 
equation in stratified media. In order to get such a regularity estimate we 
interpret equation div( 7 V?i) = 0, for a complex valued coefficient 7, as a 2 x 2 
differential system for real valued functions. If we denote by = 5Rw and 
u( 2 ) = 3it we have that the vector valued function (u^\ u^) : fl — )• R 2 satisfies 
the system 

"«* <« 

where we used the convention of repeated index summation and where 



"y = oShkSij - eS hk (SiiSj 2 - S[ 2 5ji), 

for l,j, h,k G 1,2 and with a = 9? 7 and e = 5J 7 . By assumption ((3]), system ^ 
satisfies the strong ellipticity condition 

A-w<4^<Aie 
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For this type of systems, Li and Nirenberg proved in |LN| a regularity result 
that we state here for our particular equation. The following proposition is a 
special case of Proposition 1.6 in |LN| . 



Proposition 3.1 Let and^ 2 ' be two complex constants satisfying |3J). Let 
r > and let h, g%, . . . , g n be complex valued functions of class C°° ■ 
Let v G i? 1 (i? r ) be a solution to 

div + (7 (2) - 7 (1) )l s +) Vu) =h + divg in B r , 

where g = (gi, . . . ,g n ). 

Then, for every multi-index (3' , D^v G C°(B r ) andv G C°°(B r ). Moreover, 
for every S > and k > 0, 



V W<k-i \p\<k J 



where k = k + [2=i] + 2 and C = C{8, k, n, A). 

Note that, by / G C°°(B^) we intend that / is separately C°°(B^) and 
C°°(Br). 

Observe that, as a consequence of this result, since v is continuous in B r and 
Vv is bounded separately in B~^ l _ 5 y and in fl,"jn r , then G L°°(B( 1 _ S } r ) 
and 



I|V«|U- (B(I _ 0P )< C \\v\\ L2{Br) + £ P&ftlUw+Ell^ 

'.9lU 2 (s r 

V l/3|<fc-l l/3|<fc 

(7) 

for jfc = 3 + . 

These regularity estimates can be extended to C 1:Q interfaces and less regular 
h and <?. (see Theorem 1.1 in |LN] ) . 

Our proof of Lipschitz stability estimates follows the approach used by 
Alessandrini and Vessella for the conductivity equation ( |AV] ) . In their proof a 
crucial role is played by the Green function for the conductivity equation with 
bounded leading coefficient. In our case, to our knowledge, the existence of a 
Green function in the whole domain Q is not known for an L°° complex coef- 
ficient in dimension n > 3. Existence of such a Green function is established 
in [DM] for the 2-dimensional case or for uniformly continuous coefficients in 
any dimension (see also HK for a generalization of such a result to unbounded 
domains) . 

A Green matrix for a strongly elliptic operator 
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is a matrix valued function G : {(#, y) : x, y G £1, x ^ y} — >• M™ xn such that 

LG{-,y) = S„Id in SI 
G(-,y) = ondft 

where Oy is the Dirac distribution concentrated at y and Id is the identity matrix 
in R™, 

If we have a Green matrix G = {Gij}fj =1 for system ((6]), it is easy to see 
that the first row of G, G = (Gn, G21) interpreted as a complex valued function, 
is a Green function for the operator div(7V-) in SI, in the sense that 



f 7 VG(-,y)V<l> = <I>( 2 ,) 
J ft 



for any complex valued function $ 6 Cg°(f2). 

Let us now state an existence result for the Green function for equation 
div^Vu) = for a continuous complex valued coefficient 7. 

Proposition 3.2 Le£ 7 G C (O) satisfy assumption (0) and ie< d^ = dist(x, <9S1). 
There exists a unique function G(x,y) continuous in {(x, y) € SI x SI : 1/ y}, 
locally integrable with respect to y per every x G SI and smc/i <nai, /or every 
f G G£°(S1) tfie function 



u(x) = / G(x,y)f(y)dy 
belongs to H l (Sl) and satisfies 

— div (7V 'it) = /, 

m i/ie weaA; sense. Moreover, 

7 VG(-,2/)V0 = /or ewer*/ G G °°(S1), 



and, /or every r\ G G^°(S7) swen £/ia£ 77 = 1 in B r (y) for some r < d y , 

(l- V )G(;y)eH 1 (n). 

Furthermore 

G(x, y) = G(y, x) for every x, y G SI (8) 

and 

G(-,y) <Cr 1 ~%, for every r < d v 12. 

Proof. This result follows from Theorem 5.4 in [HK] and the observation 
that, given a Green matrix for system (J6j) , it is possible to get existence of the 
Green function for equation div (jVu) = 0. The symmetry result (|SJ) follows 
from Theorem 1 in [DMj . □ 
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These results on the Green function have not been extended to L°° coeffi- 
cients, hence we cannot use a Green function for our problem. For this reason 
we will construct some solution of our equation that has the same behavior of a 
Green function, but only for y in a certain special subset of 0. Before doing this 
we need to extend our original domain ft to a Slo by adding an open cylinder 
D whose basis is the fiat portion Si of d£l fl D\ and with height greater than 
r . Let Kq = {x G D : dist(x, Si) > r /2}. If we set fl — fl U Do then dClo 
is Lipschitz continuous. 

We extend any complex coefficient 7 defined in by setting it equal to 1 in 
Do- For simplicity we will still denote this extension with 7. 

Let us consider any subdomain in fl and let us consider the chain of domains 
connecting it to D\ (see assumption (H3)). For simplicity let us rearrange 
the indices of subdomains so that this chain corresponds to Z?o, D\, . . . , Dm, 
M < N. Let us denote by S = l->fL Dj and consider a connected subset K C S 
with Lipschitz boundary such that K fl dDj = Ej ; U £j+i for j = 1, . . . , M, 
K D K a and dist(K,dS \ {S M +i U D }) > r /16 . 

Let us denote by T(x, y) the standard fundamental solution for the Laplace 
equation, given by 

r(x,y) = - ( \^\x-y\ 2 - n inR",n>3. 

where u n is the volume of the unit ball in n dimensions. 

Let 7, 5 e C. Then a straightforward calculation shows that 

!±r(x,y) + sT(x,y*) if x n >0,y n >0, 
(j + s^T(x,y) if x n -y n <0, 
±T{x,y)+tr{x,y*) if x n <0,y n <0, 

where y* = (y 1 , . . .,y n -i, -y n ), s = and t = , is a fundamental 

solution for the differential operator 

div((« R n +7l R ») V-) in R". (9) 

Proposition 3.3 Let 7 satisfy assumptions (H1)-(H3) in ri and let S, K 

and K be defined as above. For y e K there exists a unique function G(-,y), 
continuous infl \ {y} such that 

f 7 VG(;y)-V<l> = cP(y), V<f>eC^(n). 

Furthermore, 

\\G{;y)\\ H , {nXBr{y)) <Cr^ n l\ Vr<d y /2, (10) 

and 

G{x 1 y) — G(y,x) for every x,y £ K. (11) 
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Let Di and Di + \ be two subdomains ofS such that dDi + \C\dDi contains aflat 
portion satisfying assumption (H3). Let us fix the origin at Pi+i G 

and let v be the outer normal to Di at the origin. Let y = —rv for some r € 
(0, ro/6) and let x G B ro / e C\Di + \. There exists a constant C — C(ro, A, n, A, L) 
such that 



G(x,y) 



71+71 + 1 



71+71 + 1 



< c 



(12) 



V x T(x,y) 



< C 



Proof of Proposition \3.3[ Let us construct G by taking advantage of the 
fundamental solution for operator ([9|). If y € K, there is a couple of contiguous 
domains of S such that y G S D (Di U -D;+i) and dist(j/, Dj) > ^ for every 
j G {0, . . . , M} \{l,l + 1}. Let us fix the origin at the point Pi+i G Let 
us denote by T t (x, y) = r 71 7l+1 (a:, y) and by 7 = 7 - 7; lb 



7; + i1r" and let 



G(x,y) = Ti(x,y) +w(x,y), 



where w is solution to 



div ('-/V x w(-, y)) = div h in f2 

w(-,y) = -Ti(-,y) on <90, 

where h = ; yV x Ti(-,y). This problem has a unique solution because — L;(-,y) G 
H^idti) and since 7 = in fl,U%, div h G iT" 1 ^). 
Now, in order to get (|12l) consider 



where r;(:r, y) 



w (x) = w(x) +fi(x,y), 
(x)Ti(x, y) with <j> G C 00 ^), = in B ro/3 (y) and = 1 in 



^\£ 2ro/3 (2/). 

The function wq satisfies 



div (7 V wo) 
w 



div(/io + h) in SI 
on dn, 



(13) 



where ho = 7V x Li(-, y). Multiplying equation (TT3|) by wo, integrating by parts, 
using Schwartz inequality and the fact that 7 = in Di U -D2+1 and VL; = in 
B r /3{y) we get 



C 



< C||V™ ||^ 2 



7Vfz(-,2/) • VuJ 
1/2 

L 2 (n)' 
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where C = C(ro, n, A, A, L). By Poincare inequality this also implies that 

\\w \\m(n) < C (14) 



Estimate (|T0|) follows then immediately from (|T4|) and from the behavior of Tj. 

Since h = in Di U Z?; + i and /i belong to C°° f-Bro)> we can apply to 
function wq estimate ([7]) and get 

l|Vu;o||L^(B ro/6) < cl\\w \\ L * {Bro/3) + Pf<M^ 0/ 3)l ( 15 ) 

V l/3|<fe / 

where fc = 3 + [^]. 

By the previous inequality, (|T4")) and (fTS"]) we get that 

||Vw || L oo (Bro/6) < C and ||«J |U~(s ro/6 ) < C 

Finally, since w = w -fi and ||Vf i\\ L °°(B ro/6 ), ||f i|U~(B ro/6 ) < Cj C2) follows. 

Symmetry (|11[) of G in K follows by standard arguments based on integration 
by parts (see for example [EJ Theorem 13, p. 35] ). □ 

Solutions to our equation are harmonic in each of the subdomains Dj, hence 
they are piecewise analytic in tt. Now we want to show that is that we can 
analytically extend each analytic portion through the flat interface. This prop- 
erty of u will allow us in Proposition 3.5 to derive Hoelder estimates of unique 
continuation in K. More precisely we will prove the following 

Theorem 3.4 Let D\, . . . Dm be a chain and let u be a solution to 

div(^\7u) — in Dj U Dj + i 

for some j £ 1,2, ... , M . Then there exist two positive constants C\ = Ci(A, n) 
andC = C(A, ro, A, L, n) such thatu\ D can be extended by a function u analytic 
( c \ 

in the set Dj U Ej + \ , where 

E^\ ] =(xe D j+1 : dist(x,Bm{x ) n £j+i) < ^ 



INL~(^u^>) ^ CMviDiUE&y ( 16 ) 
The proof of this Theorem is contained in the Appendix. 

Proposition 3.5 Let K and Kq as before, and let v 6 H [K) be a solution to 

div (7V1;) =0 in K, 
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such that 

\\v\\l°°(k ) < C £o r 2 -", 

and 

\v(x)\ < C(e + Eo)r 2 dist(x, Sju+i) 1- * for every x € K. 

Then 

Hx)\<C Ur^— ) (E +s )r 2 - n - , (17) 

w/iere z = P M +i - 2rv {P M +i), r E (0, ±r ), r = Si € (0, 1), JVi and 



C depend on ro,L, A,n, X, and t t = — \ 3 3 ' 1 _ 2 Y , /or r\ — ro/4Ci and Ci as in 

H^r-J 

Proof. By Theorem [331 the function v\ Dn can be extended analytically to a 



function vq on Do U E[ 1 such that Vq\ d = v\ D and 



IKII £ oo (A , UB co 1 ) ) <C(£b + eo )rg- n . 

Let r\ = j^r where C\ and e[ Ci ' are the same as in Theorem 13.41 Note that, 
C\ » 4 because A > 1. 

Let us consider the sphere B^ ri of radius 4ri = ff- < ^ strictly contained 
in Kq. Let B 3ri and B ri be spheres concentric to B iri and of radius 3ri and r± 
respectively. Let P 6 Ej such that dist(P, Pi) < Let us construct a chain 
of spheres of radius r\ such that the first is B ri , all the spheres are externally 
tangent and the last one is centered at P — 2r\V\ where V\ is the exterior normal 
vector to Si. We choose this chain so that the spheres of radius 4ri concentric 
with those of the chain are contained in Dq U E± . Such a chain has a finite 
number of spheres that is certainly smaller than N\ = + 1. 

By the three sphere inequality (see, for example |ADB| ) we have that 

| Uo |l . . x^r" 1 ^ , _ M-r N i„2-n 



>{B 3ri (P-2r lVl )) ^ C£ '{Eo+Eq) 1 



where C = C(A, r ,n). 

In particular, for every P 6 Si such that dist(P, Pi) < ^P-, we have 

H«olU~(B ri (P)) < Cef 1 (E + e q ) 1 - tNi rl~ n . 
Since Sftwg and Qvq are harmonic, 

'|V«blL_y ^Cr^ef'iEo+so) 1 -^' 1 
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This also implies that the Cauchy data of v on Si are small, 

< Cef 1 (E + e ) 



E^-nfl' (Pi) 



l-r"i2-n 
'0 ' 



< Cel \E +e Q Y- r Vo 1 "". 



E-ns',. (Pi) 

V 2 

Let us now consider v\ n . Due to the transmission conditions 



dv 



l~> E+ns; ■ (Pi 



1-t"i 2-n 
'0 i 



By Trytten Cauchy estimates for solutions to elliptic equations (cfr. [T] and 
[ABRV] ) we get that 



/ 

2 16 



\Vv\ 2 dx < Cel T % * (E + ff0 ) 2(1 - rJVl5l) rf " n) 



(18) 



for some <5i G (0, 1) depending only on the a priori data. 

By standard regularity estimates for harmonic functions we get 



l~ B 1 ns m _ m (Pi 

2 8 



Again we can apply the three sphere inequality considering r\ < S with P4 ri C 
£>i and P. ri C Di fl Bra._m and applying again Theorem 13.41 we get that the 
analytic extension v\ of v\ D satisfies 



INI 

|V«i| 



L°°(s n (P)) ^ Ge o (^o + eoj 



'o > 



L~(B ri (P)) 



(E + e ) 



l-r ai iiJ-fl 



for every P G S 2 such that dist(P, P2) < ^. 

Hence again by smallness of Cauchy data, transmission conditions and Tryt- 
ten inequality ([18]) . we get 



\v(x)\ < Ceq 



(M + l)N 1 fiM + 1 



i_ t m« 1(5 « 2 _ 



for every x such that x = P — 2tr\v{PM+i) where P G Sm+1i dist(P, Pm+i) < 



^ and 1< t < 2. 



Now, let x = Pm+i — tv {Pm+i) where r < r 1 . We can again use a three 
spheres inequality for the spheres centered at Pm+i — 3r 1 v(P M+ i) and of radii 
ri , 3ri — 2r and 3r x — r and get (fTT|) . □ 
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4 Proof of Theorem 12.11 

Let -Do, D\, . . . , Dm be the chain of domains such that 

Il7 (1) " 7 (2) IIl~(^) = E := || 7 (1) - 7 (2) IIl~(o). 

Consider S, K and Kq as defined in the previous section. 

For y € K, let Gi(x,y) and G-2(x, y) be the singular function related to 
7*- 1 -* and 7^, respectively, whose existence and behavior has been shown in 
Proposition 13.31 

Let U = Q and U k , for k = 1, . . . , M, be given by U k = fi \ \J k j=1 Dj. Let 
also denote by W k = U^ =0 Dj. Define, for y, z G K, 

S k (y,z);= f {^ 1 \x)- 1 { - 2 \x))V x G 1 {x,y)-V x G 2 {x,z)dx 
Ju k 

Observe that, by (ITUl) . 

\S k (y,z)\ < CE{d(y)d(z)) 1 ^^ , for every y,z€K, 

where C depends on the a priori assumptions and d(y) = dist(y, U k ). 
The main point of the proof consists in showing that 

Proposition 4.1 For every y,z £ K we have Sk(-,z),Sk{y,-) G Hl oc (K) and 
div(j^VS k (-,z)^ =0, cfe (V^VSfcO, •)) = misT. (19) 

For the proof of Proposition ^ . 1 1 we need the following approximation result that 
we will prove later on. 

Lemma 4.2 Let 7 satisfy assumptions (H1)-(H3) and let Dq, Qq, U k and K 

as above, and, for y £ K , let G(x, y) be the function defined in Proposition \3.3[ 
Assume that {7ft,}^ 6N is a sequence of complex valued continuous coefficients 
that converges to 7 in L s (S}q) for every s e [l,+oo) and let Gh{x,y) be the 
Green's function for Then 

lim sup\\G h (;y)-G(;y)\\ m(Uk) = 0, (20) 

for every K CC K . 

Proof of Proposition [^TTj We consider, for j = 1,2 a regularization of 7^) 
obtained by convolution with mollifiers ah G C 3O (R") such that < an < 1, 
ah = for \x\ > i, J Rn ah = 0. Then consider 

lh\ x )=i l b) {v) a h{x-y)dy. 
Clearly 7^ is smooth and 9^7^ > t an d l7fc I < ^- 
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Let Gj,h be the Green's function of div ryjf^V-J. 
Let K CC if. By ([20]) 



sup||G Vi (-,y) -Gi(-,2/)|| ffl([/fc) -» as h -> +00. (21) 

Set 

S k . h {y,z) = f (j^ - 1 W)\7 x G hh (x,y)\7 x G 2 (x,z)dx, 
Ju k 

for y, z £ K. Let z £ if; by using the symmetry of G^ft and differentiating 
under the integral sign, we have 

div (7ft V5fc,/ l (-, z)) = in if, for every h. 

By US]), 

Sfc,h(-,«)->5 fc (.,«) inL°°(if). (22) 



Let ifi be such that if C if 1 C if 1 C if. By Caccioppoli inequality 16.11 we 
have that 

\\SkM{-i z )\\m(k) ^ C\\Sk,h(-,z)\\ L 2( Kl) . 
On the other hand, it is easy to check that 

||Sfc,fcMI|LW < G(dist(if 1 ,if)d(z)) 1 -^ . 

This implies that, by considering a subsequence, Sk,h(-,z) converges weakly in 
H 1 ^) to Sk(-,z) (by In particular this implies that 



div ( 7 (1) VS'fe(-,z)) =0 in ifj. 



Since ifi is any domain compactly contained in if, and since we can proceed in 
the same way with respect to z, we can conclude that (fTTJ| holds. 

Proof of Lemma \4-2\ Since 7ft —> 7 in L s (D,q) for every s > 1, there exists a 
subsequence, that we continue to denote by {7ft} that converges a.e. to 7 in fi - 

For some / & C°°(^Iq), let uq and Uh in £/q(^o) be solutions in fi to 
div(7Vuo) = — / and div(7ftVwft) = — / respectively. Observe that 

div (7ftV(«ft - u )) = div ((7 - 7ft)Vw ) in fi . 

Multiplying the above equation by Uh — uq and integrating by parts we get 

7fc|V(uft - u )\ 2 = / (7 - 7ft)Vu • V(itft - w ). 



Using the strong ellipticity condition and Holder's inequality we get 

/ \V(u h -u )\ 2 <\( f | 7 -7ft| 2 |V Mo | 2 V (|VK- Uo )| 2 ) 5 
Jn \Jn / 
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hence 

||V(uh-uo)||^ (no ) <aQT |7-7/ l | 2 |Vu | 2 ) 

By the dominated convergence theorem, we get 

\\V(u h - wo)||l 2 (o ) as h -> 
and, by Poincare inequality 

- wollff 1 ^) as ^ 0. 
Multiplying equation div(7/ l Vw/ l ) = — / by Uh we get 



n„ 



7h.|VUft| = / /lift. 

By strong ellipticity and Schwartz inequality we get 

|Vu h | 2 < A]|/]|i,2 ( n )|KI|i3 (no) 
and, by Poincare inequality 

\^u h \ 2 < C\\f\\ L 2 {no) \\\7u h \\ L 2 {no) . 
Hence 

Klliri(no) < C|l/IU-(n„) for /* > 0. 

Functions uo and Uh satisfy in K the assumptions of Proposition [3~T1 hence, 
for h > 0, 

sup K|, |V« fc | < C (\\u h \\„ m + ll/H c2+[ ^] (no) ) < G ll/ll c2+ [^] (no) 

This implies that, for some subsequence of {uh} 

Uh — > uq uniformly in K (23) 

and uq is continuous in K. 

By the properties of functions G and Gh we have that 



u o(y) = G(z,y)f(z)dz for y E K, 
J n 

and 

u h{y) = G h (z,y)f(z)dz for y E fl . 
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uniformly with respect to y G if, 

G h {z,y)f(z)dz -4 / G(z,y)f(z)dz 



By 



Let Q be a smooth domain such that Uk CC Q C ilo such that dist(Q, if) > 0. 
Now, let x G Q; by symmetry of G^ and G, 

div^VG^x,-)) =0 in if, 
div( 7 VG(x,-)) =0 in if. 

Again by Proposition 13.11 G and Gfc, satisfy a G 0,1 (if) bound uniformly with 
respect to x G Q and ZieN. 
Hence 

||G4-,2/)-G(-,y)|| i2(0) 

satisfies a G 0,1 (if ) bound uniformly with respect to ft. G N. 
Thus, there exists a sequence {y/i} C if such that 

sup \\G h (-,y) - G(-,y)\\ L2(Q) = \\G h (; y h ) - G(-, y fc )llz>«3) (24) 
ye it 

and y/j — >• y G if (as a matter of fact this holds for some subsequence to which 
we restrict). Therefore 



\\G h (;y h ) - G ft (-,y)|| L2( Q ) as ft + 



DC. 



(25) 



By ([Ml) and (gSJ) we have that, for every / G G°°(f2o) such that supp(/) C Q, 
we have 



(G h (z,y h )-G(z,y))f(z)dz 



< 



(G h (z, y h ) - G h {z, y)) f(z) dz 



(Gh{z, y) - G(z, y)) f(z) dz 



— > as ft — >• +oo. 



(26) 



As a solution of div (j h \7Gh{- 7 yh)) = in fi \ if , by [C] Theorem 9.1, p. 95] 

\\G h {;y h )\\ w u v(Q) <C (27) 

for some p > 2 and G independent of ft. Hence by Sobolev imbedding theorem, 
taking into account (|2l)|) and ([2"T)l . up to subsequences, 

G h (-,y h ) ^ G(-,y) m L 2 (Q). 

Now observe that, for ft — > +oo 

sup ||G h (-,y) -G(-,2/)|| l2( q ) = y ft ) - G(-, j/jOH^q) < 

ye/? 

< \\G h (;y) - G(;y)\\mQ) + \\G(;y) ~ wOIIl*(<3) -> 0. 
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Finally, by Caccioppoli inequality 

lim sup \\G h (-,y) - G(-,y)\\ H u Q) = 0. 

By the uniqueness of G this holds for every subsequence of the original sequence 
Gh so the claim follows. □ 



Proposition 4.3 If for a positive Eq and for some k € 1, . . . , M we have 

\S k (y,z)\ <r 2 - n e for every (y,z)GK xK , (28) 

then 

where x = Pk+i + ru^, < r < i>k outer unit normal to dDk at Pk+i and 



Proof of Proposition ^ .3\ Let us fix z £ Kq and consider v(y) = Sk(y,z). By 
Proposition ^. II we know that v solves in K equation div (7 1 Vu) = 0. Moreover, 
by (HUD 

Hv)\ < CEr'S dfetfc/, Efc+i) 1 ^ for y e W k , 

and by (gSJ 

\v(y)\<r 2 - n e foryeKo. 
Thus, we can apply Proposition 13.51 to v getting for < r < ^Hi ; 

a 1 / „ \ (i-i)(i-Tv) 



1^)1 <(^J (^o)r^-j ^ (29) 

where we denoted by fi = r( fe+1 ^ A ' 1 (5f +1 r r , Now let us consider 

v(z) — Sk(x, z) for z e If. 

which is solution of 

div (V 2) Vu) =0 in if (30) 

and satisfies 

|«(*)| < CE (rdist(z, E fc+1 )) 1- ^ for z e K. (31) 
By Proposition 13.51 and taking into account (|2U)) , ([3T)|) and (|3ip we have 
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□ 

Proof of Theorem \2.1\ Let us denote by e := ||Ai — ^W^h 1 / 2 ,h- 1 / 2 )- Let 

h ■= Il7 (1) - 7 (2) L-iw h ) for * = 0, 1, • • • , M. 

Note that, for y, z £ Dq, due to the extension to the complex case of Alessan- 
drini's identity (see formula ((33]) in the Appendix) 

(7 W -7 (2) )VG 1 (-, 2 /)VG 2 (-,z) =< (A 1 -A 2 )G 1 (;y),G 2 (;z) > 

n 



and, for i/,zG Ko, Proposition 13.31 yields 

\S k ^(y,z)\<Crl- n (s + S k ^) (32) 

Let Pk £ Sfc and y r = z r = P^ + rt/ when z/ is the outer normal vector to 
and r £ (0, ro/2). Let us write 

S k -i(y r ,y r ) = h + h, (33) 

where 

/i = / ( 7 (1) -7 (2) )VG 1 (-, 2 / r )VG 2 (-, 2 / r ) 



h = / ( 7 (1) -7 (2) )VG 1 (-,y r )VG 2 (-, 2/r ). 



From Proposition 13.31 we have that 

\h\ <Cr 2 - n E 
where po — On the other hand, again by Proposition 

I'll = l7 (1) -7 (2) l / |VGi(-,»r)||VG 3 (-,Vr)| 



Now by (I3"2l and by Proposition 14. 31 we derive that, for < r < 
\S k -i(y r ,y r )\ < C LIXT+e ) (e + 5 k -i+E)rl- n ( L-\ (34) 



Hence by (|3"3 > (|3"3 )l we get 
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and 

| 7 i 1) -7f ) |<C(E + 4-i+^) 



x \ {r kN ^8lr r ) 2 , \ n-2 



e + 5k-i+EJ \r 



-1/4 



By taking r — \\n y e ^ k ! yqrg ) J we have 

tf»-7fl.C, e + ^ 1+£ )(,„(^±^)) . (35, 

Let 

|l OK t|-(n-2)/4 for Q <t< J_ 



w(*) = 



(n-2)/4 for i > J_ 



This function is increasing, concave, lim^o uj(t) = and the function x — > 
xuj{1/x) is increasing. Inequality (1351) can be written as 

5 k < C(s + 5 k -x + E)u ( i±j*=i- V (36) 

Notice that the above choice of r is possible only if fin ^ e ^ k h ~+ E ) ) < 
but, if this is not the case, inequality is obviously satisfied. 
Since <5o = iterating (|5f)|) we obtain 

^ + £<(C+lf(£ + eK 



e + EJ ' 

where w/j is the composition of a; with itself k times. Now we recall that E = 5m 
and, hence, 

E + e < (C + l) M (e + E)u>m ( 

\E + E 

where C = C(A, A, M, L). Now, either E < e and this proves Lipschitz stability, 
or E > e and we can write 

E<2(C+l) M Eu JM (^), 



from which 



and 



that is 



< ui M (—^ 



2(C+1) M ~ J " \2EJ ' 



" !,L U(C+1) M ) " 2E' 



E < 



2 W M ^ 2{C+1) M ) 

Finally observing that 2u;^ f 1 ( 2(0+1)^ ) — -^jy 1 ( 2(0+1)" ) ^ ne c l a i m follows. 
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5 Concluding remarks 



Remark 5.1 Observe that in our result we can replace the full Dirichlet to 
Neumann map with the local Dirichlet to Neumann map. More precisely, let 
E be an open portion of dtt containing a flat open subset. Let iJc/ 2 (E) = 
{<ft G H x ' 2 (d£i) : supp (ft C £} and define the local Dirichlet to Neumann map 

1 /2 

in the following way: for (ft £ Hc£ (E) let 

< A?<ft,ip >= / jVuVv forfrip GH^ 2 (E) 
Jn 

where u solves equation div(~f\7u) = and u — (ft on dfl and v G H (Q.) such 
that v = ip on dCl. 

We observe that in our proof we apply the Dirichlet to Neumann map to 
functions whose support is contained in a neighborhood of the flat portion of 
dfl. Hence in Theorem \2.1\ we can substitute {5)) by 

- 7 2 ||l- ( o) < C||A? - Af[| £(JJ v2 )if -v 2) , 

where 

II^IUc^V- 1 ' 2 ) = su p{< A?0,^>: (f>,^eH^ 2 , \\<ft\\ H i/2 = 1, U\\ H m = l} 

Remark 5.2 We expect that Lipschitz continuous dependence of the admittivi- 
ties on the data still holds replacing the flatness condition on the interfaces with 
C ,a regularity. In fact the key ingredients in our proof are the construction of 
a singular function G{., y) for y G K , the unique continuation estimates and the 
C 0,1 estimates in K obtained by an application of the results contained in [LNJ. 

In the case of C ' a interfaces it is possible, proceeding similarly to what 
done in [AV], to make a C 1,a change of variables which straightens locally the 
interface and to prove the existence of the Green's function in the new variables. 
Also unique continuation estimates can be obtained in this case by means of three 
sphere inequality and finally C a estimates in K and local C 1,a estimates can be 
derived using the results in [LN]. 

6 Appendix 

6.1 Caccioppoli inequality 

For reader's convenience we state here Caccioppoli result for admittance equa- 
tion. The proof is standard. 

Proposition 6.1 Let u be a solution for 

dw(7Vtt)=0 in CI, (37) 
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and let Br(xo) C VL. There exists a constant C depending only on X such that 

C 



|V«| z < n> , , 

B p (x ) \K-P) JBr(x ) 



(38) 



for every p < R. 



6.2 Proof of Theorem 15731 

Let us fix xq G Sj+i such that dist(xo, Pj+i) < ^ and consider Bfj(xo) with 

i? = m 
-11 4 . 

By Proposition (J57TJ) (and the more general result in [LNj ) . -D^,m is contin- 
uous in Dj and Dj+i and belongs to C x {Dj U Sj+i) and C 1 (I?j + i U Sj+i)- 

Observe that, since equation (|3T|) has constant coefficients each ff, u G 
C(Br(xo)) and it is a solution to the same equation for any multi-index /3, 
hence u G C°°(B^(x )) and tt G C°°(B fl (a;o))- 

We can apply Caccioppoli inequality (|38| to Df/ u and get 



B„„{x ) 



V(i#«) 



< 



C 



(P'-P") 2 •/*,,(*„) 



for < p" < p' < R and C = C(A) of Proposition O 

Let iVo be the maximum order of derivative that we want to estimate in 
Br(x ) and define pk = R — jj^, so that p fc _i - p k = k = 0, . . . , 2N - 1- 
We have that 



ivttr < c 



that is 



B P1 (x ) 



E, 

fe = l"' B Pl( a; 0) 



21% 
R 



B PO (xo) 



du 



dx k 



< C 



2No 
R 



B P0 (x ) 



Since is solution to equation (l37l) for every fc = l,...,n — 1, Caccioppoli 
inequality holds and, consequently 

2 



-B P2 (tf ) 



C 



< C 



i? 



'pi^o) 



du 



,2jV V 



r" 



9a; fc 

2 



b po ( x o) 



Let now ov. 



2|/3|=m /j 



-Bp™ (^o) 

a m < I C 



|Z)f, u\ 2 and assume that 



27Vq 
R 



B P0 (x ) 



22 



We proceed by induction and consider 

n— 1 p 

2 = V" V" / 



a 2 



^ = E / i^'"i 2 = E E / _ 

|/3|=m+l JB p m +i(*o) \P\=m fc=l (a:o) 

< fp ^ )2 E / i^i 2 = ^(^) 2 E / i^i 2 



Ml / I |2 

11 ' It . 



U / I J B PO (xo) 

Hence, we have proved in particular that 



No 

2 



1/81=-? 

Finally, applying Proposition l3.il we have that 



M • (39) 

-BrOo) 



\\D^u\\ l ^ {Br(xo)) < CR-^\\D^u\\ L , (B } 

4 Q 1 u ■' 



Hence, by and recalling that iVo = we get 



1/2 [£J 
2\ "2" 



:(a;o) 

Observing that \f3f < 7# l /3!e l/31 one gets 



\\D$u\\l-(b 9 w)<MI3\(j±) (40) 



where M = C ^ / B (xo) \u\ 2 j and G x = en8X 2 . 
In particular, from f|40[) we derive 

101 



\\D^u(;0)\\ L ^ iB ^ Xo)) <Mf3\ i^j 



(41) 



which implies analyticity of $tu{x' , 0) and 3fit(x', 0). 

Let <j)(x') — 3?u(x',0) and set J7(x) =: 3?u(x) — <j>{x'). Then {/ satishes 

Af7 = -A<f>:=F in B^(x ) 
C/(x',0) = on B' R (x ). 

Observe now that, by standard regularity estimates (cfr. for example |GT1 
Corollary 2.36]) we have that 

\\vu\\ l ^ Br(X0)) < c (||f/||^ (fl - (xo)) + \\f\\ l ^ Br{X0)) 

8 \ 4 4 
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Analogously, the function W := Z)f, U, solves the problem 



AW = -A(Df,V) in B R (x ) 
W = on B^(x ), 



so that also 
IIVWI 



l°°{b- (x )) 



<C IIW 



lL°°(S- (x„)) 



Applying estimates (HOI) and (HTj) we then get that 



Il«(b- (so)) 

"3" 



< — rM/3! -f 



L~(B' H (x )) 



I? 4 



R 



Hence, also, 



1/3 1 



(42) 



which implies analyticity of ^|-5Ru(a;', 0). A similar estimate can be proved for 
3m. By the fact that Km and 3m are harmonic in B + and in _B~ separately, and 
from (grj) and (SI} we have that, 



|3ftu(a;)| = 



c 



<±_ M V 

- i?4 I R 



\x - Xq\ 
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which is convergent for 



\x — xq\ = R = 



R 

2Cl 



Hence 5?m can be extended analytically in a neighborhood of xq and the same 
can be proved for 3m. Repeating the same argument for all points x € 
such that dist(x, Pj+i ) < ^ and choosing R = ^ we have proved that Km can 

be extended analytically to the set Ej + \' . 

The same holds true for 3m, hence ud } can be extended analytically to 

Dj U and the extension u satisfies ([T 



6.3 A generalization of Alessandrini's identity 

Alessandrini's identity holds for solutions to conductivity equation with real 
valued coefficients. For sake of completeness we show here that is can be gen- 
eralized to the case of complex valued coefficient. 
Let Mi and M2 be the solutions to 



div ( 7 (fc) VM fe 



in O, 



(43) 
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for k = 1,2 respectively and let us consider the Dirichlct-to-Neumann maps 
A 7 (k) that, from now on, we will denote by A&. 
From (|43l) we have 

/ 7 (1) Vmi • VTJ=< A lUl ,v >, VveH 1 ^). 

Jn 

If we take v — U2 , we derive 

/ 7^^Vmi • Vit2 =< A1U1, U2 >, 

and, analogously, 

/ 7^ 2 - l Vu2 ■ Vui =< A 2 u2,ui > . 
Jn 

Let us show that < AiU\,U2 >=< AiU2,Ui >. Let w be solution to 

f div(7( 1 )Vw) = in Q 
\ w = U2 on 

then, 

< AiU2,ui >= / 7' 1 'Vw • Vui =< KiU\,U2 > . 
Jn 



Hence 



/ ( 7 (i) _ 7 (2))Vui • Vu 2 =< (Ai - A 2 ) u 2 ,U! > . (44) 
Jn 
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